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Substructure Synthesis: A Controls Approach 

Peiman G. Maghami * and Kyong B. Lim 1 * 
NASA Langley Research Center, Hampton, VA 23681 

Substructure synthesis from a controls perspective is considered. An efficient and 
computationally robust method for synthesis of substructures is developed. The method 
considers the interface forces/moments as control inputs and redefines the synthesis 
problem in terms of obtaining a constant gain compensator which would ensure the con- 
nectivity requirements of the combined structure. Orthogonal similarity transformations 
are used to provide simplified synthesized dynamics of the combined system. The sim- 
plicity of form in transformed coordinates are exploited to effectively deal with modeling 
parametric and non-parametric uncertainties at substructure level. Uncertainty models 
of reasonable size and complexity are derived for the combined structure from those in 
the substructure models. 


Introduction 

Substructure synthesis is a well established area in 
the modeling of flexible structures as is evident in the 
voluminous literature of which a sample of tutorial 
papers are given in [1-6]. It is concerned with the 
modeling of the structural dynamics of substructures 
(or components) and then synthesizing them to pre- 
dict the combined structural response. This synthesis 
is accomplished by enforcing deflection compatibility 
and force equilibrium at all substructure interfaces. 
There are several advantages of modeling via sub- 
structure synthesis, among these, (1) it allows much 
independence in the design and analysis of substruc- 
tures, which is especially helpful if substructures are 
designed, fabricated, and even tested by different or- 
ganizations. An example is the damping synthesis for 
the Space Shuttle 7 and more recently the component 
modules for the International Space Station which are 
built by different companies in several countries. (2) 
it increases the power of existing finite-element anal- 
ysis and design programs by allowing analysis/design 
by components especially in problems where too many 
finite element degrees of freedom are required to per- 
form a dynamic analysis/design of the complete sys- 
tem. (3) It allows a direct synthesis of substructure 
test data. This is particularly useful for very large 
structural systems (such as the International Space 
Station) that cannot be tested as a whole. This also 
means that it can be used as a part of an experimental 
verification tool for substructures before deployment 
as a connected structure. 

Since its earliest work in, 1 the primary goal of 
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substructure synthesis methods has not changed, i.e. 
to accurately predict the combined modal parame- 
ters: structural resonant frequencies and mode shapes. 
Consistent with the primary goal, most attention has 
been given to the issue of selecting a most effective sub- 
set of component (or assumed) modes and the related 
issue of modeling the substructure interfaces. Thus, 
relatively little attention has been given to issues ger- 
mane to active control. 

In this paper, we re-examine the substructure syn- 
thesis method from a controls perspective. In particu- 
lar, our end goal is to develop a physically realistic and 
highly efficient way of modeling and controlling large 
flexible structures based on substructure (subsystems) 
models. As in any studies in multivariable control, 
we are concerned with the issues of closed loop sta- 
bility and performance robustness due to inevitable 
modeling errors, configuration changes, or exogenous 
disturbances. Since the stability and robustness is 
strongly dependent on the degree of accuracy of the 
mathematical model to the physical structure, we also 
examine the issues of modeling uncertainties for the 
substructures and their influence on the output of the 
connected structure. 

Substructure Synthesis 

For simplicity of presentation, consider two compo- 
nents, shown in Fig. 1, that have to be synthesized. 
Consider the interface forces and moments (between 
the two components) as control input forces yet to 
be determined. Moreover, assume that collocated 
and compatible with the interface forces /moments 
are linear/angular velocity outputs from each com- 
ponent. With these considerations, the linear and 
time-invariant dynamics of the two components are 
represented as follows 

comvonent 1 
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Actual Structure 



Substructure Configuration 


Fig. 1 Synthesized substructures 


M\Z\ + D\Z\ + B\Z\ — Biu -f II \ i>\ (1) 

Vi — C ri z\ ( 2 ) 

Upi — ^d\ Z\ -f- L ri z\ 4- L ai z\ (3) 

component 2 

M 2 Z 2 + D-2Z-2 + K-2Z2 = -B 2 U + H 2 p 2 (4) 
y 2 =C r2 Z 2 (5) 

y p 2 = Ld 2 z 2 + L r2 Z 2 + L a2 Z 2 (6) 


where, for component i , Zj, denotes the displacement 
vector; M t is the positive definite inertia matrix; 
D,j is the damping matrix; Kj is the non-negative 
definite stiffness matrix; Bj is the influence matrix 
for the interface forces/moments; u is the interface 
force/moment vector; y t is the interface velocity vec- 
tor; C ri is the corresponding output influence matrix; 
y Pi is the performance output vector; which could be 
a combination of displacement, velocity, and acceler- 
ation outputs; and L^, L ri , and L ai are the corre- 
sponding influence matrices. Since the order of a large 
flexible space structure (LFSS) can be quite large, for 
design and analysis purposes the order of the system 
is reduced to a design size using model reduction tech- 
niques such as modal truncation or modal cost analysis 
to obtain for component i, i=l,2 

M ri q ri + D ri q r . -f K ri q ri = ±4>, 7 Bpa + <t> 7 //,/;, 

( 7 ) 


y Pi — ^di^iQri "b ^ri *&iQn (9) 

where q ri is the vector of modal amplitudes ob- 
tained from transformation Zi = ; M u , D n , 

K ri are, respectively, the generalized inertia, damp- 
ing and stiffness matrices; and <T>, is a matrix 
whose columns are the r open-loop eigenvectors as- 
sociated with the included modes. If the mode 
shapes are normalized with respect to the inertia ma- 
trix, and modal damping is assumed, then M ri = 
I r x r i Dn = Diag{ 2(i i uji i , . . . ,2( ri ui ri }, and K n = 
Diag{u>l }, where ojj. and Qa j = 1, • • • ,r 

are the open-loop frequencies and damping ratios for 
the ith component. Defining the state vector Xj = 
{ Qn Qn } T tl = 1,2, the second-order dynamics of 
the components can be rewritten in first-order forms 
as follows 

component 1 


x\ = A\X\ + Biu + Hipi ( 10 ) 

i/i = Cixi (11) 

y Pl = L 1 x\ + D Ul u + D pi pi ( 12 ) 

component 2, 

x 2 = A 2 x 2 + B 2 u + H 2 p 2 ( 13 ) 

y -2 = C 2 x 2 ( 14 ) 

y P 2 = L 2 x 2 + D U2 u + D P2 p 2 ( 15 ) 


where A t , B t , H,, Ci, and Lj, i = 1,2, respectively 
denote the state matrix, input influence (interface 
forces or moments) matrix, exogenous disturbances in- 
fluence matrix, velocity output influence matrix, and 
non-interface output influence matrix for each com- 
ponent. The matrices D Ul , D Pl , D U2 , and D P2 , 
represent feedthrough matrices associated with non- 
interface outputs of each component; xi and x 2 denote 
the state vectors of the components; y\ and y 2 denote 
the interface velocity output vectors; y Pl and y P2 de- 
note the non-interface output vectors components; and 
pi and p 2 represent the exogenous disturbances acting 
on the components. The system matrices in the first- 
order forms are related to those in second-order forms 
as follows: 


V i — C ri 4 * ,q r; 


( 8 ) 


Aj = 


0 I 
l\ , —D r 


Bi = 


±<f>i T Bi 


(16) 
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Hi = 


$; T Hi 


c ri = [ o r\«i>, ; ; 


(17) 


The dynamics of the system in new coordinates be- 
comes 


Bi [ -Ldi*&i L ai Q %K ri L Vi ^i L ai <$>iD ri J 

(18) 



B Ui +/„, < l> l < l>/ / ih ; /!,,, ' //, (19) 

Now, for the two components to be connected, in a 
dynamical sense, one has to find a control input vector 
u, representing the internal forces and moments at the 
interface, such that the displacements and rotations of 
the two components at the interface remain identically 
the same for all p\ and p 2 , and all compatible initial 
conditions. To achieve this, append the dynamics of 
the components, to obtain 



y = y i - y-2 = [ -c 2 ] ( } (21) 


where p = 

A\ A-2 

-4.3 ,4.4 


{£} 


and 


= [ n c r c y 


Ai 0 

0 A 2 


[ N c R c ] 
(26) 


51 

5 2 


= [ N c R c ]- 


Si 

So 


(27) 


Si 

So 


= [ So: Rc y 


Si 0 

0 So 


(28) 


Note that the linear and angular velocities at the in- 
terface are collocated and compatible with the forces 
and moments, i.e., C\ = Sf and — Co = Bj , such that 
along with Eq. (23), one has 


Up — 


[ Vrn \ _ 

' Li 

0 

f ;Ci 

} + 

' D U1 " 

u+ 

n t 

' ' 

\ y P 2 J 

0 

1/2 

\ X-2 

J 

Du,2 


c 

B 2 


= Sf[Ci -Co ] T = 0 


(29) 


D 


Pi 

0 


D 


V 2 


{£} 


(22) resulting in 


The new output vector y is the difference between the 
velocity of the two components at their interface. The 
problem is to find u such that y is identically zero for 
all pi and p 2 , and compatible ;ri(0) and .t 2 (0). From 
Eq. (21), for y to remain identically zero for all p\ 

and p- 2 , j ' | must remain identically in the right 

null space of the matrix [ Ci — C 2 ] . This means 
that an appropriate control vector u must be found 
that would render the system identically unobservable, 

| identically remains in the undetectable 

subspace of the system given in Eqs. (20) and (21) for 
all pi and p 2 . Let N c denote an orthonormal basis for 
the right null space of [ C i — C 2 ] , i.e., 

[ Ci -C 2 ] N c = 0 (23) 

and let R c denote an orthonormal complement to N c . 
Transform the system in Eq. (20) via an orthogonal 
similarity transformation, such that 

{ Z } = [ ^ i } < 24 > 



Bi = JVf 


Bi 

B 2 


= 0 


(30) 


Now, the control input vector u must be chosen such 
that f3(t ), which represents the coordinates incom- 
patible with the connectivity of the two components, 
remains identically zero for all t, p, and ’’compatible” 
initial conditions, i.e., /3(0) = 0. To accomplish this, 
the input vector u must be chosen to render j3 unreach- 
able from p, that is 

u = -B^[A 3 a + H 2 p\ (31) 


Note that the feedback gain associated with a and 
the feedthrough term associated with p must be in the 
form given in Eq. (31) to render (3 unreachable from p. 
Although, any u in the form of u = — B T 1 [A 3 a+H- 2 p} + 
R/3, with R as a non-destabilizing arbitrary matrix, 
would also be feasible, the condition for synthesizing 
the dynamics of the combined system necessitates that 
3 = 0 for all t, thereby reducing the expression for u 
to that given in Eq. (31). The matrix inversion in 
Eq. (31) is guaranteed, as long as both input influence 
matrices Bi and ll> are full rank, which they are since 
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the interface forces are distinct. Therefore, since B\ 

By 


and B -2 are full rank, so is 


So 


, and so is B 2 , since 


B-2 = B'r 


By 

b 2 


, and R c is an orthonormal basis for 


the column space of 


By 

b 2 


. Using the control input 


u of Eq. (31) in Eq. (25), with the aid of Eq. (30) 
simplifies to 


r a 1 

r Ai 

Ao 1 

f a 

\ 

' Hy ' 

11 

0 

l ^ 

l p . 

r 

0 


(32) 


It is obvious from this equation that if the initial con- 
ditions are compatible, such that /3(0) = 0, then f3 
remains identically zero for all p and t, which means 
that, in the physical coordinates, y\ = y 2 for all p 
and t. Therefore, compatibility conditions between 
the two components would be satisfied for all p and 
t. The synthesized dynamics of the combined sys- 
tem can be written in a reduced-order form in the 
transformed coordinates, or in the original component- 
based coordinates. In the transformed coordinates, 
these dynamics are given as follows 

a = Aid: + HyP (33) 


with the non-interface output defined from Eqs. (12), 
(15), (24), and (31) 


Vp = { V * ) = Nia + D lP 
l yp 2 J 

where 


Ah 


Dy 


' Ly 

O 

/V. - 

' D U1 ' 

0 

To 


c 

Du,2 

f B Pl 

0 



D U1 ' 

O 

Dp 2 


Du,2 


B^A 3 


B^H 2 


(34) 


(35) 


(36) 


The dynamics of the combined system may also be 
presented in the original coordinates. Redefine the in- 
terface feedback law (Eq. (31)) in terms of the original 
coordinates, to obtain 


Using Eq. (37) in Eqs. (20) and (22), results in the 
combined dynamics of the two components in the orig- 
inal coordinates. 


{£} 


Ay-ByGy -ByG 2 
—B 2 Gy -4 2 — 5 2 G 2 

Hy — ByDy —ByD 2 

— S0D1 H 2 — B 2 D 2 


]{:} 
]{::} <s 


+ 


yp — 


By — D Ui Gy —D Ui G2 
—D U2 Gy 


D Pl D Ul D y 
—D U2 Dy 


l 2 — D U2 G 2 
— D Ui D2 

D P2 Bu2 B 2 


{::}+ 
{*:} <39) 


Model Uncertainty 

In this section, issues of parametric and nonpara- 
metric uncertainties in the component models and the 
way they effect the synthesized dynamics of the system 
will be considered. Parametric uncertainties include 
uncertainties in modal frequencies, damping ratios, 
and mode shapes of the components. The nonpara- 
metric uncertainties considered here are the unmod- 
eled dynamics of the components. If modal models are 
used to represent the component dynamics, nonpara- 
metric uncertainties would be due to the truncated 
modes of the components. The development in this 
section uses the representation of the synthesized dy- 
namics of the system in the transformed coordinates, 
as given in Eq. (33). This would make the treatment 
of parametric and nonparametric uncertainties easier. 
Here, only the case wherein there are no feedthrough 
terms (acceleration terms) associated with the non- 
interface outputs (see Eq. (34)), is considered. The 
treatment for the general case would be similar but 
more involved. 


Uncertainty in Frequency and Damping 

As mentioned earlier, the uncertainties in modal fre- 
quency and damping are embedded in the elements of 
the state matrices A\ and . One can represent these 
uncertainties via additive uncertainty blocks AAi and 
A Ao , such that the true state matrices are given by 


u = -B^[A 3 Nj I ^ } + H. lP ] 

s -[o, Gs ]{*;j- [Dl ^]{«} 

(37) 

with the constant gain matrices Gy and Go, and the 
feedthrough terms Dy and D 2 defined as 

[Gy Go ] = B^A 3 Nj ■ [Dy D 2 ] = B^ 1 H 2 


A 1 -a A 1 + AAi (40) 

Ao — > Ao -f- AAo (41) 

Note that no assumption is made on the structure of 
the uncertainty blocks, A A\ and AAo, i.e., they can 
be structured or unstructured. Using Eqs. (40)-(41) 
in Eq. (33), and assuming that there are no uncer- 
tainties in the mode shape data and nonparametric 
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uncertainties present, one obtains the uncertainty rep- 
resentation in the transformed coordinates. 


Up = {Ni + LM c A)a 


= (Ai + A^4 j)q; + Hip 


where 


where 


AAi = N / 


N T C1 AAi N C1 + Nl A A 2 N C2 (43) 

A r Cl and N C2 are column partition of N c according to 
the components. This indicates that the additive para- 
metric uncertainties in the component dynamics in the 
original coordinates carries through to the synthesized 
system in transformed coordinates, in the sense that it 
still remains in the form of an additive parametric un- 
certainty. However, it is observed from Eq. (43) that 
the structure of the uncertainty does not remain the 
same, its size decreases to match Ai, but its bound re- 
mains invariant, since columns of N c are orthonormal. 
In other words, AAi would have the same bounds as 
that of A A, i.e., 

\AA\<6 -A |A4i| <6 (44) 

The parametric uncertainty problem given in Eq. (43) 
may easily be put in a linear fractional representation 
(LFT) form for robust control design and/or analysis. 
Also, one could observe from Eq. (38) that, in this 
case, treating the uncertainty in the transformed co- 
ordinates would be cleaner and easier than treating it 
in the original coordinates. 

Uncertainty in Mode Shapes 

The parametric uncertainties in the mode shape 
data affect the elements of the input and output in- 
fluence matrices, which include matrices C i and Co. 
This means that there would be uncertainty associ- 
ated with the elements of the basis vectors N c , as seen 
form Eq. (23). Here we assume that the uncertainty in 
the mode shapes, whether structured or unstructured, 
can be translated into an uncertainty representation 
for N c . With this assumption, the actual basis vec- 
tors N c may be represented in various representations 
of uncertainty, such as additive, multiplicative, etc. 
Here, let N c be represented as follows 

N c -A N c + M c A (45) 

where M c is a chosen basis for uncertainty influence, 
and A represent the uncertainty, which can be struc- 
tured or unstructured. Using Eq. (45) in Eq. (33), 
and assuming that there are no uncertainties in the 
frequency and damping data or nonparametric uncer- 
tainties present, one obtains 

a = (ii + AAi)a + (Hi + A T Mj H)p (46) 


AAi = A ' M ' A M. A + A T MjAN c + Nj AM c A 


Hi 0 


Hi 0 

0 H-2 


Li 0 

0 To 

(49) 


From these equations, it is observed that the input and 
output parametric uncertainties in the dynamics of the 
components, in the original coordinates, diffuse into 
input, output, and state parametric uncertainties in 
the transformed coordinates representation of the syn- 
thesized system dynamics. Moreover, the state para- 
metric uncertainties involve a quadratic form of the 
uncertainty A. This form of uncertainty can be dealt 
with the technique presented in [8-9], wherein poly- 
nomial uncertainty representations are represented in 
optimal linear fractional representation (LFT) form 
approximations. The remaining uncertainties given in 
Eqs. (46) and (47) may easily be put in a linear frac- 
tional representation (LFT) form for robust control 
design and/or analysis. 

Unmodeled Dynamics 

The traditional approach in substructure synthesis 
has been to use higher bandwidth in the dynamics 
of the components to be synthesized, roughly twice 
the bandwidth of interest for the combined structure. 
This approach has worked successfully in many appli- 
cations. However, it is somewhat ad hoc, with some 
potential shortcomings. First, it does not provide 
guaranteed levels of accuracy for the predicted system 
parameters, such as frequencies, damping ratios, etc. 
Second, it may suffer from finite element method’s po- 
tential loss of numerical accuracy for higher frequency 
modes in addition to a loss of correlation with mea- 
sured frequencies. As mentioned earlier, control design 
requires accurate assessment of model parameters or 
potential uncertainties. Traditionally, unmodeled dy- 
namics, which typically include the higher frequency 
modes not included in the design model, have been 
treated via additive uncertainties in the plant model, 
which forces the control system to roll off to avoid 
potentially destabilizing spillover problems. Unfortu- 
nately, unmodeled (truncated) dynamics at the com- 
ponent level generally do not correspond to the same 
in the synthesized (combined) system. Consequently, 
it is imperative that the effects of unmodeled dynamics 
in the components be characterized at the system level 
for proper dynamics and controls design and analysis. 
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The approach taken here is similar to the treatment 
for parametric uncertainty, which essentially consid- 
ered the effects of uncertainty in the transformed co- 
ordinates. Rewrite the dynamics of the components, 
including the unmodeled dynamics, in a first-order 
form. 

component %= 1 or 2 


[ X/ 1 

' Ai 

0 

f Xi 


' Hi ' 


' Hi ' 

U r 

0 

~Ai 

u 

r 

. B > . 

u + 

. H i . 





(50) 

Vi — [ C'i 

Ci ]{ 

Xi 1 
Xi J 

(51) 

Upi = [ Li 

L, ]| 

C Xi | 

l I 

(52) 


Here, the overbars indicate terms associated with the 
unmodeled dynamics. It should be noted that in most 
applications an accurate knowledge of the parameters 
associated with the unmodeled dynamics is not avail- 
able, thus they are typically characterized in the form 
of uncertainty. Also, note that, similar to the previous 
treatments, the feedthrough terms in the non-interface 
outputs have been omitted. Following the synthesis 
procedure outlined previously, the component dynam- 
ics, as given in Eqs. (50)-(52), are appended to obtain 
a system similar to the one given by Eqs. (20)-(22), 
except that the order of the states are rearranged such 
that the modeled dynamics’ states appear first, fol- 
lowed by the unmodeled dynamics’ states, i.e. , the 
state vector is given by [ xf x\ xj xZ ] ■ Then, 
the space of feasible states (in the connectivity and 
compatibility sense) are characterized by an equation 
similar to Eq. (23), but for the expanded system, as 
follows 


where 


Ai = N? 


Hi = N, 


Ai 

0 

0 

0 

An 

0 

0 

0 

Ai 

0 

0 

0 


' H x 

0 

n t 

0 

It- 

c 

Hi 

0 


L 0 

H- 


0 

0 

_0 

An 


Nr 


(56) 


(57) 


Ni = 


L\ 

0 


0 

Ln 


Li 

0 


_0 

Ln 


N r 


(58) 


It is reasonable to expect that both [ C i —Cn ] and 
[ Ci —Cn Ci —Cn ] are full rank. Define the 
null space of the matrix [ Ci — Cn ] by the matrix 
N Cl , whose columns are orthogonal. Then, it can eas- 


ily be shown that the columns of matrix 


N C1 

0 


is 


included in the null space of [ Ci —Cn C i —Cn ], 
represented by the matrix N c , from Eq. (53). Now, 
choose and partition N c , such that 


N c 


N C1 Nc j 
0 N- C2 


(59) 


Here, the first column partition corresponds to the 
states of the modeled dynamics, and the second parti- 
tion to the states of the unmodeled dynamics. More- 
over, the first row partition corresponds to the null 
space of matrix [ Ci —Cn ] . Using Eq. (59) in Eqs. 
(54) and (55), and separating the states, gives 


— -^1 -A-2Qt u H - H m p 

(60) 

&u — ^3 CX-m + ^4 O^u + H u p 

(61) 

Up — NmOt-m “1“ N u Ot u 

(62) 


[ Cl -Cn Cl -Cn ] N c = 0 (53) 

Following the coordinate transformation (similar to 
Eq. (24)) and deriving the interface input vector u 
to guarantee the compatibility of the components’ dis- 
placements and rotations at the interface, results in the 
synthesized dynamics of the combined system, which 
in transformed coordinates takes the form 


where 


A i 0 


Mi 




i 

CO 



O 

(M 


L o o 

N C1 N- Cl 
0 N- C2 


0 0 
_0 0 
Ai 0 

0 A 2 

(63) 


a = Aid + Hip 


y p = Nia 


(54) 

(55) 


' H m ' 



N- Cl ' 

. H u _ 


L 0 

N- C2 _ 


Hi 

0 

0 

H- 

Hi 

0 

0 

77 - 
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Fig. 2 Synthesized cantilevered beam problem 


[ N m N u ] 


Li 

0 

Li 

O 


£ 

0 

L-2 

0 

L-2 \ 


' CO 

(M 

£ 
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Note from Eq. (60) that A\ is the same as A\ of the 
reduced dynamics for the baseline combined model 
(see Eq. (33)), so that if no unmodeled dynamics 
were present, Eq. (60) would represent the assem- 
bled structure. Furthermore, the state equations for 
a m and a u are coupled. One way of looking at this is 
to imagine that there are two systems, one attributed 
to the modeled dynamics (a m states) and the other 
to the unmodeled dynamics (a u states), which are in 
feedback connection. As mentioned earlier there are 
uncertainties associated with the unmodeled dynam- 
ics parameter, namely, Ao, A3, A4, H u , and N u , i.e., 
they are not known accurately. It can be said that 
the uncertainties, due to unmodeled dynamics, ap- 
pear in feedback connection to the nominal plant. In 
other words, the uncertainties associated with unmod- 
eled dynamics in the components, which are typically- 
represented by additive uncertainty to the component 
models, takes the form of an uncertainty in a feedback 
loop with the nominal combined structure, when the 
synthesized dynamics are considered. 


Numerical Results 

The proposed methodology for substructure synthe- 
sis and controls is applied to a two-dimensional can- 
tilevered beam problem to assess its effectiveness and 
feasibility. The current application involves substruc- 
ture synthesis and does not include the substructure 
uncertainty treatment. The cantilevered beam is as- 
sumed to consist of two substructures, a cantilevered 
beam and a free-free beam, as indicated in Figure 2. 
The two interface forces/moments as well as the ver- 
tical displacement and slope are indicated, along with 
an exogenous disturbance at the tip of the second sub- 
structure. A performance output representing the tip 
deflection is also considered. The material and geo- 


metric properties of the system were chosen to provide 
considerable modal content in the low-mid frequency- 
range to make the synthesis and control design task 
more challenging, and are provided in Table 1. The 
structure and its two substructures were modeled as 
Euler-Bernoulli beams. Each of the two substructures 
was modeled using the first 40 flexible modes, resulting 
in an 80th-order state space model for each of the sub- 
structures. The 40 modes used were deemed sufficient 
to provide a basis for describing the first 20 modes 
of the combined structure. The combined structure 
was also modeled using its first 20 modes, resulting in 
40th-order state space model. This model was used 
for validation of the synthesis approach. 


Table 1. Geometric and structural properties 


Property 

Substr. 1 

Substr. 2 

Whole 

Length 

150 

50 

200 

Mass/Length 

1000 

1000 

1000 

Rigidity, El 

1.0e7 

1.0e7 

1.0e7 

No. of modes 

40 

40 

20 


Using the procedure outlined in Eqs. (23)-(39), the 
dynamics of the two substructures were synthesized 
(combined), resulting in a 160th-order model in the 
original coordinates or a 158th-order model in trans- 
formed coordinates. Note that the synthesized model 
may be further reduced in order by using model re- 
duction techniques, such as modal cost analysis or 
Hankel norm techniques. Table 2 compares the natu- 
ral frequencies of the synthesized model vs. that of the 
validation (truth) model. Comparison of the frequen- 
cies indicate a good match between the two models, 
thus demonstrating the feasibility- of the proposed syn- 
thesis technique. It should be pointed out that the 
degree of matching between any two frequencies in the 
table depends on whether the modal content used in 
the dynamics of each of the components provides a 
sufficient basis to express a mode at the system level. 
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Table 2. Natural frequencies 


Mode No. 

True 

Synthesized 

1 

0.0088 

0.0088 

2 

0.0551 

0.0553 

3 

0.1542 

0.1559 

4 

0.3022 

0.3053 

5 

0.4996 

0.5004 

6 

0.7464 

0.7477 

7 

1.0425 

1.0522 

8 

1.3879 

1.4006 

9 

1.7827 

1.7850 

10 

2.2268 

2.2315 

11 

2.7203 

2.7469 

12 

3.2631 

3.2939 

13 

3.8553 

3.8600 

14 

4.4968 

4.5075 

15 

5.1877 

5.2407 

16 

5.9279 

5.9856 

17 

6.7175 

6.7252 

18 

7.5564 

7.5769 

19 

8.4447 

8.5348 

20 

9.3823 

9.4767 


In order to obtain a better assessment of the appro- 
priateness of the synthesized model, the Bode magni- 
tude plots of the transfer function from the exogenous 
disturbance to the tip displacement were computed 
and are shown in Figure 3 for the synthesized and truth 
models. The two transfer functions agree very closely, 
as observed from Figure 3. Note that the peaks of the 
magnitude plots are finite although undamped mod- 
els were used to compute the transfer function. The 
reason for the contradiction is the coarseness of the 
frequency points, i.e. , no frequency point was chosen 
exactly at a resonance. 

These results demonstrate the feasibility of the pro- 
posed synthesis approach. 

Concluding Remarks 

A novel and robust method for synthesizing the 
dynamics of substructures has been developed. The 
interface forces and moments were considered as con- 
trol inputs, and used to design a constant gain com- 
pensator which synthesizes the combined dynamics of 
the substructures. The synthesized structure can be 
realized in original coordinates, or alternatively, in 
a transformed coordinates (obtained via orthogonal 
transformations). The realization of the system dy- 
namics in transformed coordinates has the advantages 



Fig. 3 Bode magnitude plot from exogenous dis- 
turbances to tip displacement; true(solid), synthe- 
sized (dashed) 

of being more amenable to uncertainty and robust- 
ness analysis and design as well as being smaller in 
order. Using this realization, uncertainty models of 
reasonable size and complexity have been derived for 
the combined structure from parametric and nonpara- 
metric uncertainties in the substructure models. These 
models can and will be used for robust control design. 
An application of the proposed synthesis technique to 
a cantilevered beam problem demonstrated its feasi- 
bility. 
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